Two routes to magnetic order by disorder in underdoped cuprates 
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We study disorder-induced magnetism within the Gutzwiller approximation applied to the t — J 
model relevant for cuprate superconductors. In particular, we show how disorder generates magnetic 
phases by inducing local droplets of antiferromagnetic order which eventually merge, and form a 
quasi-long range ordered state in the underdoped regime. We identify two distinct disorder-induced 
magnetic phases of this type depending on the strength of the scatterers. For weak potential 
scatterers used to model dopant disorder, charge reorganization may push local regions in-between 
the impurities across the magnetic phase boundary, whereas for strong scatterers used to model 
substitutional ions, a local static magnetic moment is formed around each impurity. We calculate 
the density of states and find a remarkably universal low-energy behavior largely independent of both 
disorder and magnetization. However, the magnetic regions are characterized by larger (reduced) 
superconducting gap (coherence peaks) and a sub-gap kink in the density of states. 

PACS numbers: 74.20.-z, 74.25. Jb, 74.50.+r, 74.72.-h 



I. INTRODUCTION 

Local experimental probes have highlighted the impor- 
tance of intrinsic disorder and spatial inhomogeneity in 
the studies of high-Tc superconductors. In particular, 
scanning tunneling spectroscopy (STS) measurements 
have revealed nano-scale modulations in the gap for 
quasiparticle excitations and the local density of states 
(LDOS)>i _ — Complementary to these (energy resolved) 
density modulations, it is well-known that the spin den- 
sity is also heterogeneous and exhibits spin-glass behav- 
ior in the underdoped regime. This seems to apply to 
both "clean" cuprates like YBa2Cu306+a: (YBCO) where 
quasi-static SDW order is found in the far-undcrdoped 
regime^— and to intrinsically disordered materials like 
I^-zSr^CuC^ (LSCO) where the static spin correla- 
tions are long-range, and persist for a large doping range 
well into the superconducting domeJ^— The size of the 
spin-glass phase in temperature and doping is clearly en- 
hanced by disorder. At present a detailed understanding 
is lacking of the connection between the modulated spin 
density and the STS tunneling dataJ^ 

In order to further probe the effects of disorder, 
samples with additional extrinsic impurities have also 
been studied. Most notably, substitutional Zn ions 
strongly suppress T c and induce localized low-energy 
density of states resonances^ In the spin sector, 
Zn is known to slow-down and eventually freeze the 
spin fluctuations i 13 ' 1 7 ~ — For example, in near-optimally 
doped LSCO increasing amounts of Zn substitution has 
been shown to shift spectral weight into the spin gap, and 
eventually generate elastic magnetic peaks in the neutron 
response^ A similar Zn-induced spin-freezing has been 
obtained for YBC O 20 ' 21 and is generally consistent with 
^SR data on underdoped cuprates j^ 3 - - — 

Theoretical modeling of these experimental results are 
complicated by the need to include both charges, spins, 



and realistic disorder configurations j22r— Within an un- 
restricted Hartrce-Fock approximation of the Hubbard 
model including <i-wave pairing one may acquire a quali- 
tative understanding of the disorder-induced magnetism; 
due to the splitting of low-energy in-gap impurity states, 
it can be advantageous to generate local AF order which 
may percolate, and eventually form a quasi-long-range 
ordered state . 31 i 35 ' 36 This scenario is a many-impurity 
generalization of the well-known one-impurity result of 
induced local magnetization; 3 ^— and agrees with trans- 



port measurements 
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Recently, the dynamics of spin 



freezing by impurities, i.e. the shift of magnetic spec- 
tral weight to low energies, were calculated explicitly and 
shown to be consistent with this scenario as well4£ 

The above approach, while successful to an extent, has 
obvious theoretical drawbacks. First, it cannot describe 
the approach to the Mott insulator; underdoping has to 
be understood as the effective increase of correlations rep- 
resented by U/t as one underdopes due to the suppression 
of screening, but there is no way to calculate this effect 
systematically. Secondly, the connection between corre- 
lations and pairing is entirely artificial, since the BCS 
pairing term is added by hand and treated in mean field. 
What is needed is a technique which allows for the study 
of local variations of observable quantities in the presence 
of disorder, which easily accounts at least for the crude 
effects of strong correlations in the underdoped regime. 

Here, we study disorder-induced magnetism within 
the Gutzwiller approximation (GA) of the t — J model. 
The GA has not been traditionally applied to include 
spin non-degeneracy, but was extended to include anti- 
ferromagnetism by Ogata and Himeda for the homoge- 
neous case . 47 ' 48 The so-called extended Gutzwiller factors 
(EGFs) have been applied to inhomogeneous problems 
by interpreting them in a site-dependent manner. The 
site-dependent EGFs have been used to study local mag- 
netic moments around a nonmagnetic impurity^ 7 - elec- 
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tronic states around a vortex core^ and anti-phase su- 
perconducting domain structures^ More recently, a sim- 
plified version of the EGFs, with the advantage that they 
reduce to the well-defined extensively tested original site- 
dependent Gutz wilier factors^ have been used to exam- 
ine the energetics of the charge and spin stripe ordered 
superconducting stated We focus on two distinct cases: 
1) weak scatterers with impurity concentration equal to 
the doping level, ni mp = S, modeling the disorder poten- 
tial from out-of-plane dopants, and 2) small concentra- 
tions, n imp = 1 — 2%, of strong scatterers simulating the 
effects of substitutional disorder e.g. Zn or vacancies. In 
both cases, we hnd that disorder may induce AF phases 
in the undcrdoped regime. The origin of the induced 
magnetism is vastly different, however. In the former 
case of dopant disorder, the redistributed charge density 
creates local regions close to half-filling in-between the 
impurities, pushing these areas across the bulk magnetic 
phase boundary. We denote this kind of magnetic phases 
by type I. By contrast, in the other case of substitutional 
disorder, even a single impurity nucleates magnetization 
in its vicinity, and the final spin structure consists of over- 
lapping single-impurity regions which may orient them- 
selves in order to minimize the exchange energy^! This 
kind of disorder-induced magnetic phase is called type 
II in this paper. The density of states in the disordered 
phase largely retains its d-wave 'V shape in agreement 
with experiments 4 - - — and consistent with earlier theoreti- 
cal studies i 32 i 44 ' 53 i 54 Surprisingly, however, the magnetic 
regions exhibit a larger superconducting gap and reduced 
associated coherence peaks in agreement with the general 
characteristics of the large-gap regions in the experimen- 
tal STS data. Finally we discuss additional sub-gap fea- 
tures in the LDOS which could function as fingerprints 
for local magnetism in tunneling experiments. 



II. MODEL 



The t — J Hamiltonian, defined on a 2D lattice, is given 



by 
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where c\ a creates an electron at site i with spin tr. 
is the spin operator for site i and Pq is the Gutzwiller 
projector defined by P G = JJ^l - ?V™4), where h ia = 
cycler is the spin dependent number operator. For all 
results shown in this paper we have used J/t = 0.3 and 
t' = —0.25t. In Eq. ([T]), (ij) denotes nearest-neighbor 
sites with associated hopping amplitude t whereas (ij) 
refers to both nearest- and next-nearest neighbor sites 
with hopping amplitudes t and t' , respectively. Disorder 
is introduced into the system by N point-like scatterers 
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To solve the t — J model, the no double occupancy con- 
straint has to be approximated. Zhang el alJ& intro- 
duced the Gutzwiller approximation (GA) to replace the 
Gutzwiller projectors in Eq. (JT|); in this paper we use a 
simplified version of the EGFs giving rise to the following 
rcnormalizcd Hamiltonian 
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The simplified EGFs depend on the local values of the 
magnetic and pairing order parameters, the local kinetic 
energy, and hole density defined by 
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where |\&o) denotes the unprojected ground state wave 
function. The simplified EGFs are given as 
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Note that 

the simplified EGFs allow for A-j- ^ Aj.. The rewriting 
in Eas.([5|l-([T5]l of the EGFs is identical to that used by 
Yang et al£2- 

A direct diagonalization of the Hartree-Fock Hamilto- 
nian Hh-f obtained from a mean-field decoupling in Eq. 
([3]) is not sufficient because the simplified EGFs also de- 
pend on the order parameters. Instead, the energy has to 
be calculated from the mean field Hartree-Fock Hamilto- 
nian and then minimized with respect to the unprojected 
wave function \^o) under the constraints of both fixed to- 
tal electron density n, = N e , and fixed wavefunction 
normalization (^ol^o) = 1^ This is equivalent to mini- 
mizing the function: 



W = (* \H H - F \*o) - A«*o|*o> - 1) 




(14) 



3 



which leads to the following renormalized mean-field 
Hamiltonian 

H mi= £ 7h~ (15) 



with the self-consistent equations 
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Here tr denotes the opposite spin of a. The derivatives 
of the EGFs entering these equations can be straight- 
forwardly derived from Eqs.(|g])-( [T3"]) . We have solved 
these unrestricted equations self-consistently by iter- 
ation on 24 x 24 lattices, by diagonalization of the 
Bogoliubov-de Gcnnes (BdG) equations associated with 
the excitation operators 7^ and 7 ncr defined by 6V|- = 

E„ (M7nt + < it i)^ 

At low doping, such a simple iterative procedure fails 

. t 

to converge in the inhomogeneous case because -3-^- 
diverges when the local density approaches half-filling, 
<5 = 0. Therefore, small changes in the electron den- 
sity between consecutive iterations causes large jumps in 

, 9i , In this paper, we therefore restrict the doping level 
of the inhomogeneous cases to be larger that S > 0.115 
where iterations are sufficient for obtaining numerical 
convergence. On the other hand, for the homogeneous 
case allowing for a superconducting and a commensurate 
AF phase, the phase diagram can be easily mapped out 
for 5 > 0.026 and is shown in Fig. [TJ T his pha se diagram 
is similar to that obtained e.g. in Rcfs. I48ll57t Note that 
the phase boundary for the magnetic order in Fig. [1] is 
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FIG. 1: (Color online) Phase diagram for the homogeneous 
case with commensurate (tt, tt) AF order m, d-wave supercon- 
ductivity A, and kinetic energy \ [ see Eqs. (|4l6|) ]. 



valid only for a standard (tt, tt) AF order. At finite dop- 
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FIG. 2: (Color online) Type I behavior. Magnetization shown in a real-space field of view for doping levels close to the 
magnetic phase transition of the clean system 5 = 0.115, 0.125, 0.13, 0.135 (left to right). The black dots show the positions of 
the impurities. The dopant disorder is modeled with a potential strength Vi = t, and n; mj) = 8 where ni mp is the impurity 
concentration and 5 the doping level. As seen, the dopant-induced magnetic regions gradually disappear as the doping increases. 
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FIG. 3: (Color online) Type II behavior. Magnetization shown in a real-space field of view in the presence of 1% (a,c) and 2% 
(b,d) strong impurities with Vi = 1004 for doping S = 0.14 (a,b) and 5 = 0.15 (c,d). 



ing striped magnetic order can also be stabilized within 
a similar approach^ exhibiting a slightly different phase 
boundary. 

III. RESULTS 

We begin the results section by discussing the case of 
dopant disorder, i.e. n,; mp = S where each scatterer is rel- 
atively weak, Vi = t. In previous studies of unrestricted 
Hartree-Fock applied to the Hubbard model it was found 
that each dopant induced local magnetization, leading to 
a scenario where the amount of disorder-induced mag- 
netization is proportional to the doping level contrary 
to experiments.— Such an approach, however, does not 
include any band-widening with increased doping, and 
can be made consistent with the lack of magnetization 



in the overdoped regime only by requiring U/t to be a 
decreasing function of doping. The present model nat- 
urally includes the effects of strong correlations in the 
underdoped regime and, as seen from Fig. [2] the dopant 
disorder indeed induces a finite magnetization but only at 
low doping levels. The magnetic phase shown in Fig. [2j 
which we denote type I, is incommensurate as seen from 
Fig. |4ja), and disorder- induced as verified by a vanish- 
ing magnetization in the absence of disorder (not shown). 
The origin of the magnetization is a charge redistribution 
caused by the weak impurities: it is energetically favor- 
able for the electrons to be located away from the disorder 
sites which then push these local regions across the mag- 
netic phase boundary similar to the homogeneous case 
shown in Fig. [T] Clearly, this mechanism is dominant 
at low doping where regions more readily reach the local 
critical doping level. 



The origin of the magnetization in Fig. [2] is different 
from the disorder-induced AF studied previously, which 
we denote by type II, where each impurity gives rise to a 
local magnctizationi 31 ' 37 ~ — In the single impurity case, 
each defect gradually "freezes" magnetic fluctuations as 
correlations or the impurity potential increase, shifting 
spectral weight from high to low energies, eventually ac- 



cumulating in a local u = peak. Within the present 
model, we reproduce the latter mechanism in the limit of 
strong scatterers as shown in Fig. 03 Specifically, Fig. [3] 
shows a situation with Vi — lOOt for n imp — l%(a,c) 
and Tii mp = 2%(b,d), and doping 8 — 0.14(a,b) and 
5 = 0.15 (c,d). Though care must be taken when ex- 
trapolating single-impurity to many-impurity effects in 
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d-wave superconductors ; 36 ' ' 58 ~ — the results shown in 
Fig. [3] can be explained from overlapping single-impurity 
magnetizations. The result agrees with the general no- 
tion that enough strong scatterers at a given doping level 
can lead to static magnetic order. For higher doping 
levels larger concentrations of impurities are needed to 
freeze the spins as shown explicitly by comparing e.g. 
Figs. EKa) and EJc) ^o^d 6 

Next we focus on the result shown in Fig. [2ja). Figure 
2] shows other relevant physical quantities for this param- 
eter set: (a) the Fourier transform of the magnetization, 
(b) the absolute value of the magnetization (c) the 
electronic charge density rii, (d) the gapmap extracted 
from the LDOS, (e) the peak-height extracted from the 
LDOS, and (f) the gapmap extracted from the LDOS 
without the possibility for magnetic order. A compari- 
son of Figs. [4jb,c) verifies the direct correlation between 
the local density and the induced magnetization for type 
I behavior. Hence, the origin of the magnetization is lo- 
cal phase transitions caused by charge modulations with 
regions closer to half-filling pushed across the magnetic 
phase boundary [see Fig. [1] . 

The LDOS Ni(u>) can be obtained from 

Ni(w) = Y,9ia [|tw| 2 *(w - E n ) + \v nia \ 2 5{u + E n )] , 

n 

(19) 

where i denotes the site index, n is the index of the eigen- 
states with BdG eigenvalue E n , and gf G = is the 

Gutzwillcr renormalization factor originating from the 
"hopping" between different times entering Eq. (|19|l . 53 
Surprisingly, the local superconducting gap shown in Fig. 
BJd) extracted from the LDOS is larger in the magnetic 
regions of the system contrary to the expectation from a 
conventional competitive scenario. The larger gap results 
from a combined effect of 1) lower effective doping in the 
magnetic regions and hence a larger local pairing order 
parameter [see Fig. Q] caused by the impurity induced 
redistribution of the electron-density, which is enhanced 
by the magnetic order, and 2) a renormalization of the 
local pairing constant in the presence of magnetization 
dominated by the last term proportional to mirrij in Eq. 
(jTTJ) (see also Fig. [7]below). Comparing the gapmap with 
and without magnetic order in Fig. HJd) and Fig. @|f) 
clearly illustrates that magnetic order contributes signif- 
icantly to the spatial modulation of the gap measured by 
STS in the underdoped regime^— However the two bulk 
arguments above are not sufficient to explain the entire 
modulation in Fig. HJd) indicating a nontrivial effect of 
the spatial inhomogeneity. 

The peak height of the coherence peaks also varies in 
real space and tends to be anti-corrclated with the lo- 
cal gap magnitude as seen from Fig. Ste). This anti- 
correlation is a direct consequence of g™ in Eq. (TT9")) 
which reduces the low-energy spectral weight in large gap 
regions due to their closer proximity to half-filling. The 
anti-correlation between local gap and coherence peak 
height is consistent with STM measurements^ which 
motivated a picture of locally modulated pairing in the 
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FIG. 4: (Color online) (a) Fourier transform of the magneti- 
zation averaged over 10 distinct impurity configurations, (b) 
Absolute value of the magnetization |m«|. The white dots 
show the positions of the impurities, (c) Electronic charge 
density m. (d) Gapmap extracted from the LDOS (half the 
distance to the positive coherence peak), (e) Peak height 
extracted from the LDOS (positive coherence peak), (f) 
Gapmap extracted from the LDOS, but where the possibility 
for magnetic order is removed by hand. The parameters are 
identical to those used in Fig. [2jb). 



cupratc superconductors^ - — Within the present strong- 
coupling approach, we find a similar anti-correlation for 
type 1 disorder-induced magnetization in the underdoped 
regime. 

Representative curves for N(u) in different local gap 
regions are shown in Fig. [Ha,b) where one clearly sees 
the anti-correlation between local gap and coherence 
peak height. As discussed previously, the presence of 
a robust 'V shaped density of states at low energies 
is novel and not contained within models that ignore 
electronic correlations ! 32 ' 44 ' 53 ' 54 This can be seen explic- 
itly from Fig. [5jc,d) where we compare the spatially 
averaged LDOS within the present model and a con- 
ventional Bogoliubov-de Gennes approach^ The latter 
model clearly piles up states inside the gap. By contrast, 
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FIG. 6: (a) the local chemical potential Hi in the x direction 
for a 24 x 24 system containing a single impurity with Vi = t 
(dashed) or Vi = 1004 (solid). The impurity is situated at 
site=13 and S = 12.5. (b) the same as (a) but cut off right 
before the impurity site to highlight the spatial dependence 
away from the impurity site. 
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FIG. 5: (Color online) LDOS N(uj) vs. energy u. (a) LDOS 
for two sites in Fig. H£b) (red curve is from a large-gap region 
[site (13,21)], blue in a small-gap region [site (13,11)]. (b) 
LDOS for the same two sites but in Fig. HJd). (c) Spa- 
tially averaged LDOS corresponding to the four panels in 
Fig. [2] (black 5 = 0.115, green 6 = 0.125, blue 5 = 0.13, 
and red S = 0.135). (d) The spatially averaged LDOS for 
a dirty rf-wave superconductor with t' = — 0.25t, J = Lit, 
and n imp = 0%, 1%, 5%, 10%, 20% (bottom to top at lu = 0). 
The calculations were done using 20 x 20 or 24 x 24 lattices 
with 10 x 10 supercells and (d) was averaged over 10 different 
impurity configurations. An artificial broadening of the delta- 
functions in Eq. (|19[) with r\ = 0.016 was used to smoothen 
the curves slightly. 

the spatially averaged LDOS shown in Fig. E^b) corre- 
sponding to the four panels in Fig. [2] displays a remark- 
ably universal low-energy density of states. The magne- 
tization at low doping causes additional minor structure 
inside the gap. We interpret these sub-gap kinks as a 
disordered version of a similar in-gap kink caused by al- 
tered contours of constant energies in a model with bulk 
coexistence of AF and d-wave superconducting orderi 15 i 65 
Such features could serve as tunneling fingerprints of lo- 
cal magnetism in the underdoped regime. 

Next, to underline the importance of minimizing the 
energy and to understand the difference between a strong 
and weak scatterer we define a " local chemical potential" 
Hi defined from 

dW 

^- = -fa-Vl)-IM. (20) 

Figure [5] shows the local chemical potential Hi f° r differ- 
ent impurity potentials. It is clear from Fig. [6ja) that 
Hi works against the impurity potential on the impurity 
site since a larger Vi leads to a smaller However, the 
renormalization of the impurity potential by Hi has the 



largest impact on weak impurity potentials because Hi 
takes values in the range [0, 2.15]. Figure^b) illustrates 
how Hi varies for the neighboring sites depending on the 
impurity potential. For weak impurities Hi spreads out 
the impurity potential, while for large impurity poten- 
tials Hi attracts electrons to the neighboring sites. The 

dominant contribution to Hi turns out to be -j-^- which 
is related to the kinetic energy. Thus for weak impurity 
potentials, the holes on the impurity site tend to gain 
kinetic energy by spreading out as a consequence of the 
change in the effective hopping integral, while the oppo- 
site is the case for a strong impurity. To summarize, the 
strong correlations treated within the GA, affect the im- 
purity potential in two distinct ways. First, the suppres- 
sion of the electron density at the impurity site i enhances 
the effective hopping integral t e ff = ght which acts as 
a healing of the damage done by the impurity. Second, 
the electron density redistributes to minimize the energy 
because the Gutzwiller factors themselves depend on the 
electron density. The redistribution of the electron den- 
sity is enforced by the local chemical potential which also 
tends to work against the impurity potential. 



IV. CONCLUSIONS 

In this work, we have studied disorder-induced magne- 
tization within the t — J model with correlations treated 
using the Gutzwiller approximation to implement the no- 
doublc-occupancy constraint. In general, the inclusion 
of correlations strong enough to describe band narrow- 
ing and other crude features of the Mott transition lo- 
cally was seen to justify the phenomenological descrip- 
tion of underdoping in terms of a renormalizcd Hub- 
bard U jt used until now in Hartree-Fock treatments of 
disordered correlated d-wave superconductors. However, 
some unexpected subtleties were also discovered. In the 
case of dopant disorder where the impurity concentration 
equals the doping level and the individual disorder po- 
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tcntials are weak, local charge reorganization can induce 
regions with finite magnetization located away from the 
impurities. On the other hand, strong scatterers locally 
pin AF regions which may merge and eventually form a 
quasi-long-range ordered spatial structure. Remarkably, 
the LDOS at low energies remains largely disorder in- 
dependent whereas the superconducting gap (coherence 
peak height) extracted from the LDOS is increased (de- 
creased) in the magnetic regions present in the under- 
doped regime. 

The goal of this analysis is to work towards a the- 
ory incorporating disorder, superconductivity and cor- 
relations capable of describing local spectroscopies of 
cuprates across the phase diagram. Perhaps the most 
successful approach thus far, in terms of reproducing the 
many statistical observables reported by STM, has been 
the phenomenology of Nunncr et aZ.— To some extent 
this success was only possible, however, because a) cor- 
relations were neglected and doping dependence ignored; 
and b) an impurity was assumed ad hoc to add indepen- 
dent Coulomb and pairing potentials to the system. The 
present work has attempted primarily to address the de- 
ficiency represented by a). Earlier works, notably Ref. [66| 
also incorporated disorder and correlation in the super- 
conducting state, but neglected both the modulation of 
the pairing interaction by disorder and the possibility of 
local magnetism. The latter effect has been shown here 
to lead to some of the important correlations present in 
/iSR, neutron scattering, and STM; in particular the anti- 
correlation between peak height and gap amplitude has 
been exhibited in the locally ordered magnetic phases. 

On the other hand, the modulation of the pairing po- 
tential represented by the exchange constant J in the 
model renormalized by the local Gutzwiller factors is rel- 
atively weak, as is shown in Fig. [7J Thus the positive 
correlations of the gap size with the O defect position, 
identified as crucial in McElroy et al.^ will not occur 
in the present model when the correlations become weak 
enough in the overdoped phase. From Fig. [7] one can see 
that the impurities cause local reductions of J because 
the Gutzwiller factor g s is a decreasing function of the 
doping level. The present approach, therefore, does not 
include the possibility, discussed in Nunner et al.^ that 
the defect distorts the lattice locally leading to a differ- 
ent pairing interaction, as represented, e.g. by a local 
enhancement of J. Local impurity-enhancements of J is 
also obtained by explicitly including the different impu- 
rity potentials on the two sites involved in the exchange 
process^ The notion of local J-enhancement has been 
explored by Maska et al.^ Foyevstova et al.^ Johnston 
et al. and Khaliullin et al^- with model-dependent re- 
sults; all however indicate that sizeable modulations of 
the pairing interaction can indeed occur. These effects 
have been left out of the current approach, but will be 
included in our future studies of this problem. 



For completeness, we mention a second possibility 7 -^ 
to explain the defect-gap correlations at optimal to over- 
doping, that a second source of disorder is present, corre- 
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FIG. 7: Real-space plot of the effective exchange coupling 
Ji = "Ylij where j are the four nearest neighbors to site i, 

obtained from Jij = \ 7f5~~ I ^ a shown for a configuration 
of weak impurities [type I disorder] with rii mp — 5 — 0.125 
similar to Fig. [4] 



latcd with the O dopants, and not imaged in the exper- 
iment of McElroy et al£L This scenario is plausible but 
requires the existence of a second, independent source 
of disorder. In addition, it has been found in density 
functional theory structural studies that both in the O 
dopant case and that of the structural supermodula- 
tion, Cu04 half-octahedra are tilted in identical ways 
by the perturbation , 72 ' 73 and are correlated empirically 
the same way with the gap modulations , 61 ' 74 ' 75 lending 
credence to the idea that a single set of O dopants is 
the primary driver of the structural distortions and gap 
changes. We therefore believe that the combination of ef- 
fect of electronic correlations as described here, together 
with a practical description of the modulation of the pair- 
ing gap, should provide a complete description of the 
statistics of local STM observables over the whole phase 
diagram. 
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